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ABSTRACT 

Aims. We investigate potential models that could explain why multiple proto- stellar systems predominantly show single jets. During 

their formation, stars most frequently produce energetic outflows and jets. However, binary jets have only been observed in a very 

small number of systems. 

Methods. We model numerically 3D binary jets for various outflow parameters. We also model the propagation of jets from a specific 

source, namely L1551 IRS 5, known to have two jets, using recent observations as constraints for simulations with a new MHD code. 

We examine their morphology and dynamics, and produce synthetic emission maps. 

Results. We find that the two jets interfere up to the stage where one of them is almost destroyed or engulfed into the second one. 

We are able to reproduce some of the observational features of LI 551 such as the bending of the secondary jet. 

Conclusions. While the eff'ects of orbital motion are negligible over the jets dynamical timeline, their interaction has significant 

impact on their morphology. If the jets are not strictly parallel, as in most observed cases, we show that the magnetic field can help 

the collimation and refocusing of both of the two jets. 

Key words. ISM: Herbig-Haro objects - ISM: individual objects: LDN 1551 IRS 5 - ISM: jets and outflows - stars: formation 

1 ■ Introduction studies of stellar populations, however, show that large num bers 

of bi nary and multiple star systems are expected (Patience et al. 



or bi nary ana multiple star systems are expectea (.ratience et al. 
Well-collimated jets occur across a broad range of mass scales 2002j: ISimon et al.ll995l:lGhe^ et al.lll993l:S"uquennov & M~ 



}>? from young stellar objects to active galactic nuclei. In young 1991) ^^ ^^^h systems present associated disks and jets ? It has 

rN stellar objects (YSOs), jets and outflows are more easily ob- ^een argued that the disk configuration in a binary system may 

^ served at optical, infrared and UV wavelengths and the proper- either be in circumstellar disks or c ircumbinary disks dep end- 

■ ■ tiesare relatively well known (Reipurth & Bally 2001; Ray et al. i^g on the star separation distance d HartJganetal. 1 120003) . On 

|2007|). YSO jets are beheved to be launched by means of a ^^e other hand, it is well established that multiple sources can 

magnetic field which is anchored or frozen into a circumstel- ^^ the source of HH objects, e.g. T Tau, IRAS 04 325-1419, Z 

lar disk, and pinched and twisted by the disks rotation. Once the Cma, Sz 68, SR 24, S Cr A, AS 353 (Reip urth etal . 1993). Now, 

pinched-in field reaches a critical angle, fast-rotating gas from ^ multiple system may produce a single jet or outflow or a set 

the disk may be loaded onto the field lines and thus launched into ^f multiple outflows. Although there should be "many" visible 

the parent cloud in an outflow (|Blandford & Payne 1982). The binary jets, only a few binary jets from binary protostars have 

most accepted schools of thought diverge on whether the wind ^een observed. The examples known to the authors are compiled 

is launched primarily from the^^inner disk region (disk wind, ^ TablelB The frequency of existing binary jets is low compared 

|Ferreira||1997|; |Komgl & Pudritz 2000) or the X-annulus where ^ ^^e large number of protostellar binary sources. 

the young star's magnetosphere interacts with the disk (X-wind, ^. ^. , , , ^, ^ .^. . , 

/I 01 . ^ \ Uf^f^rl \ TT • -.1 .1 • . • 1 1- J . Observations have also shown that multiple or quadrupo- 

(IShu et al.ll2000D ). However in either case the let is believed to , . ^ t^ • ^, ^ 1 1 Arknox a- j j 

/ — ; — :; — ^~~T r ., 11 1 . .1 .• 1 . .1 . lar lets may occur. Reipurt h et al.l (^1993) discovered a second 

brake the rotation of the disk and is thus essential to the star r, tttti /i/i ^ ^u ^ tttti o i?^ a ^ o j^ • ^ 

r ^' •. 11 .1 . . . . •. r^ 1 now HH144 irom the same source as HHl-z. lGredel & Reipurthi 

formation process, as It allows the star to accrete up to Its final /mn/iux u ^ ^ ^i /ttttioi\^ ~Z 

^ ^ (,,19941^) observed a second flow (HH121) from the same source 

as the well-known HHl 1 1 outflow. In both cases there are large 

In all such models a smgle star with its associated circumstel- ^^gj^, ^f separation between the two jets. This can be explained 

lar disk produces a bipolar outflow or jet. Estimates derived from ^f ^^^ ^^^^^ ^^^ ^^^ coplanar (Mathieu 1994). However, obser- 
vations of apparent double jets can also be explained by other 
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Table 1. Evidence of Binary Jets and Outflows? 



Outflow (Source) 


Natal Cloud 


RA (J2000) 


Dec (J2000) 


Reference 


HH154 (L1551 IRS 5) 


Tau-Aur 


04 3134.20 


+ 18 08 04.8 


Liseau et al. (2005) 


HHl-2 - HH144 


Orion 


05 36 22.85 


-06 46 06.6 


Reipurth et al. (19^ 


HH111-HH121 


Orion 


05 51 46.07 


+02 48 30.6 


Gredel & Reiuurth Q994a) 


L723 (IRAS 19156+1906) 


Cep 


19 17 53.16 


+ 19 12 16.6 


Anglada (2004) 


HH288 (IRAS 00342+6347) 


Cep 


00 37 11.07 


+64 03 59.8 


Guethetal. (2001) 


HH377 (IRAS 23011+6126) 


CepE 


23 03 13.9 


+6142 21 


Ladd & HodaDD (1997) 


IRAS 16293-2422 


pOphE 


16 32 22.8 


-24 28 33 D 


Hirano et al. (2001) 


IRAS 20050+2720 


CygRift 


20 07 06.7 


+27 28 53 


Bachiller et al. (1995) 



found a distinct multi-polar morphology. They concluded that 
the lobes of the jets were in fact shell walls of two cavities swept 
clear by a single bipolar outflow. Angladaet al. (1991) discov- 
ered a double radio source at the centre of the outflow struc- 
ture and this led to a reappraisal of the of the double outflow 



flieory (Angladaet al. 1996; Gira rt et alJ 119971: lAnglada et al. 



1998; Hirano etal. 1998; Palacios & Eiroa 1999; Shirlev^elaL 
2002; Estalefla et al. 2003; Anglada 2004). Other example of 
multiple molecular outflows include IRAS 16293-2422, IRAS 
20050+2720. 

In order to explore these possibilities, and to investigate 
about the nature of jets from multiple systems, we perform nu- 
merical simulations of binary jets with typical physical quanti- 
ties. Our aim is to demonstrate the efl'ects of their relative sizes 
and speeds on their interaction, in order to test if the jet evolu- 
tion can lead to situations compatible with the few binary jets 
observed, and with the corresponding observed properties. 

As a starting point, we focus on jets emanated independently 
by two sources that launch the two outflows with a small angle of 
separation. A suitable test case for this scenario is the the binary 
jet HH 154, emanating from the source L1551 IRS 5. 

Numerical models of optical jets have evo lved over the 
past tw enty years, from the original models by 'Norma n et al.l 
(119811) of the Bl andford & Rees ( 1974) t win-exhaust model, to 
later models bv iHartmann & Raymondl (Il984l) : iHartigan et al.l 



(119871) : [Raymond et"an (Il988h which were 1.5D bow shock 
models which explicitly tried to model observed fea- 
tures, to axisymmetric and 3D geometries and includ- 
ing such physics as atomic radiativ e cooling (iBlondin et alJ 
'1990; Raga & Norie ga-Crespol IT998I: IStone & Harded I200QI: 
iRagaetal...200a .20011) ^ 



In our case we model a fully three-dimensional binary 
jet system using ATLAS, a new shock-capturing, multi- 
dimensional, constrained transport, adaptive-grid, direction- 
ally unsplit, higher-order Godunov astrophysical MHD code, 
ATLAS and examine its morphology and propagation dynam- 
ics. 

The paper is structured as follows: in Section [2] we describe 
the observed properties of the binary jet LI 551 IRS 5, that con- 
stitutes our test case. In Section [3] we describe the methodology, 
the set of equations used, the microphysics included, the numer- 
ical code and the initial and boundary conditions on the com- 
putational grid. In Sect. |4] we present the results of the simula- 
tion, and finally in Sect. 5 we discuss our results and the insights 
brought forward by our model. 



2. Observations of jets from the binary Protostar 
L1551 IRS 5 

The picture being built up over 25 years of observations of the 
object L1551 IRS 5 and its associated jets and outflows is that 
of a pair of YSOs each with its own associated disk and outflow, 
the whole structure em bedded in a larger disk. 

IStrom et aP (Il976h observed the near-infrared source IRS 5 
within the Lynds molecular cloud 1 551 (L1551). The source 
has an infrared emission nebulosity (lHodapplll994 l). consistent 
with the presence of a channel perpendicular to the high-density 
disk, from which the light from the cent r al star escapes and ir- 
radiates the nebulosity. Rodrigu ez et al.l (Il998h confirmed that 
L1551 IRS 5 was a binary system and show the first images 
of the circumbinary disks and the red lobes of the two jets. 
Additionally they suggested that there is a circumbinary struc- 
ture and a large-scale envelope around LI 551 IRS 5. Lay etal] 
(1994) observed a large circum binary disk of ~ 160 AU in di- 
ameter. [Rodriguez et al.l (l2QQ3h confirmed the jet binarity at ra- 
dio wavelengths, while iLiseau et"an (I2QQ5 ) provided a value for 
binary separation of 40 AU. We add that recent observational 
evidence has suggested that the source may have a small third 
co mpanion (Lim & Takakuwa2006_) . 

ICudworth & Herbigl (1 19791) discer ned two fast He rbig- 
Haro objects near the IRS 5 source. ISnell et aP (l9Su) saw 
for the first time the "remarkable double-lobed structure" in 
L1551. This was the among the fir s t mole cular bipolar out- 
flows discovered. iBieging & CohenI (Il985h ident i fied IR S 5 
as a binary source. iMoriarty-Schieven & Wanniej (Il99ll) and 
Pound & Bally ("1991") identified a second outflow from the 
same source. Mundtetal. (1991) also observed two indepen- 
dent rows of knots although these were interpreted as edges 
of a limb-brightened cav ity, and were only identified as jets by 
iFridlund & Liseaul (Il998h (hereafter FL), who measured a jet an- 
gular separation of approximately 20°. 

Recent observations carried out with the Hubble Space 
Telescope a nd SUBARU have fu rther clarified the morphology 
of the flow (iFridlund et al.ll2005l) . The optical jets extend south- 
west and disappear at approximately 1400AU from the IRS5. 
The north and south jets appear to be launched fr om th e south 
and north disks respectively. According to Lisea u et al.l (l2005h 
the northern jet is faster with a radial velocity projected in the 
plane of the sky of ~ 430 km s"^ whereas the southern jet has 
a radial velocity at most 65 km s"^ Itoh et al. (2000) observed 
a twisted morphology in the light of iron lines, and argued that 
the precession of the sources is too slow to be responsible for 
this efl'ect. Possible other mechanisms can be magnetic in na- 
ture e.g. the Lorentz forces (JFendt & Zinneckej|l998L |2QQQ|) or 
the bending may be caused by the decrease in kinetic energy al- 
lowing the ambient magnetic p ressure to be c omparable to the 
ram pressure of the jet. Indeed, .Scarrottl (11988.) observed a high 
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Table 2. Observational Data for L1551 IRS 5 outflows 



N Jet Radial Velocity 


200-250 km s'' ^ 


S Jet Radial Velocity 


65 km s-i b 


Binary Separation 


45 AU'^ 


Angle Between Jets 


20°^ 


Orbital Period 


255 years ^ 


Ambient Density 


5000 cm-3f 


Jet Density 


500 cm-3 g 


Jet Length 


1400 AU h 



" lHartiganetaD(l2000bl) 
^ iHartiganetaPfeOOObh 

' 'Liseauetal.'('2005') 

^ Fridlund & Liseau (!l998h 
" Itoh et al. ( 200QD 
^ Fridlund & LiseaS (I l998[ ) 
g Fridlund & Liseau (1991) 
^ [Fridlund & Liseaui (119981) 



3. Methodology 

Our main goal is to model binary jets. Contrary to single jet 
propagation, that can be modelled also analitically in an axisym- 
metric approximation, the propagation and interaction of binary 
jets is an intrinsically three-dimensional problem. A numerical 
study then becomes necessary, and requires the use of a fully 
three-dimensional time-dependent model with magnetic field 
and radiative cooling. Constraints on current computational ca- 
pacity require us to make a number of approximations however, 
that lead to the adoption of an ideal magneto-hydrodynamics 
(MHD) approach. In this way we are able to include the eff'ects 
of the magnetic field without losing the advantages of using a 
simple macroscopic description for a fluid particle. We addition- 
ally include the non-ideal physics of optically thin atomic radia- 
tive cooling losses. Moreover, in our scheme the ionisation frac- 
tion of hydrogen i s followed in the sim ulation . This is importan t 
as several studies ( Bacciotti & Eislofl 'el 1999; Podio et al.|[2006h 
have demonstrated that this quantity in jets decouples from the 
thermal gas conditions and follows the evolution of gas recom- 
bination. 



degree of polarization in the optical light emitted by the region 
and concluded that it coul d be explained by a toro idal field in the 
cloud around the outflows. ILucas & Rochel(ll997i) also observed 
this "peculiar pattern of alignment". The physical cause for the 
polarimetry pattern remains unexplained. 

For completeness we mention also that X-ray emission 
from the region of the head of the no rth jet was observed by 
lFavataetal.r(l2Q02h . iBallv et al.1 (l2003h performed a higher an- 
gular resolution study, finding a source of X-rays in L1551 IRS 
5 at the location of the source which they attributed to either fast 
shocks or reflect ed x-rays from IRS 5 scattered out through the 
outflow cavity. Bally et al.l ([2003) put forward a number of mod- 
els to explain the X-ray emission, including the intriguing possi- 
blity of a quasi- stationary X-ray luminous shock maintained by 
interacting colliding winds between the two protostars. 



3.1. System of Equations 

The ideal MHD equations evolve in time the four quantities (two 
vector and two scalar) which are conserved in a volume: density 
p, momentum, pu, magnetic flux density B and energy density, 
E. 



dp 

dt 



+ V • (pu) = 



pu u + {p*) I -h B B 



^(pu).V. 

^B 

^-+V-(u0B-B0u) = O 
ot 

dE 

— + V . [(£ + /7*) U - (U . B)B] + Lcooling = 



(1) 

(2) 
(3) 
(4) 



1 arcsec 



IRS 5 



^BBl. north 



NORTHERN JET 




+1 arcsec knot z 



arcsec 



SOUTHERN 
JET 



^ 



A\ 



3.5 arcsec 



4^ 



7 arcsec 




Fig. 1. iFridlund & Liseaui (Il998h provide a HST R-band image 
of the two jets from the binary protostar LI 5 51 IRS 5, located in 
the constellation Taurus. (Image courtesy of IFridlund & Liseaui 
(fl998.) ) 



where E is the total energy, arising from the sum of three 
terms, namely kinetic, internal and magnetic energy: 



E = ip|u|2 + -^ 

2 y- 1 



.i|BP 



and /?* is the total (thermal plus magnetic) pressure: 

P* = P + \\R[' 



(5) 



(6) 



Lcooiing represents the losses due to optically thin radiative 
cooling, and is described in the next paragraph. The units are 
chosen so that B absorbs a factor of 1 / V47r. The adiabatic index 
is y = 5/3 for a monatomic gas throughout the simulations. The 
equation of state is the ideal gas equation (p=nkT). 

3.2. Microphysics 

Protostellar jets are strongly cooled by radiation of collisionally 
excited lines, many of which are optically thin. It it this emission 
which render the objects visible, particularly in Ha and in the 
light of forbidden lines like the doublet of singly ionised sulphur 
at about 6700 A([SII]/1/16716, 6731). Modelling electron transi- 
tions directly, however, is neither feasible nor desirable in the 
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fluid approximation. Therefore, we represent the microphysics 
of energy lost due to radiative cooling using a cooling func- 
tion adapted from Sutherland & Dopita (1993) and depicted in 
Figure [2l which represents on a macroscopic scale the rate of 
energy loss at our temperature range for a sample (solar) set of 
abundances. [Sutherland & Dopital (Il993l) include the efl'ects of 
electron collisional ionisation, radiative and dielectronic recom- 
bination and line radiation within their cooling function. 



NEQ COOLING CURVE: ZERO FIELD, LOG(To) =6.50, SOLAR ABUNDANCES 




4.5 5 5.5 6 

Temperature, log(T) 



Fig. 2. Log normalised radiative cooling loss plotted against log 
temperature ([Sutherland & Dopital[T993l) . 



To calculate the energy loss in the gas one also needs to know 
the free electron density, which is co niputed explicitly at each 
time step and location in the fluid, as in' Falle & Ragal (|l995 l). As 
mentioned above, this is necessary because the ionisation frac- 
tion of hydrogen is not in equilibrium with the thermal properties 
of the gas, due to the fact that at the observed densities the re- 
combination times are long with respect to the crossing time of 
the visible portion o f the jet. For a thorough discussion of this 
point see Raga etlD Sl991): Bacciotti & Eislofl'el (1999). 

The cooling function shown in Figure [2] is highly non-linear. 
If we use the cooling timestep to evolve the overall jet the pro- 
cess will be extremely time-consuming. Therefore we separate 
the cooling from the dynamical evolution at each timestep, first 
calculating the dynamical change and then computing the radia- 
tive cooling using a short substep. 

3.3. The numerical code 

We computed the three-dimensional simulations using ATLAS, 
a new modular, parallel, shock-capturing, directionally unsplit, 
adaptive-grid, multi-dimensional, constrained transport, higher- 
order Godunov astrophysical MHD code. Validation and verifi- 
cation tests, including but not limited to the Orszag-Tang MHD 
vortex, the Brio-Wu and Ryu- Jones suite of shock tube tests 
and 2D and 3D blast wave tests, have been run against the 
code to build confidence in its ability to form correct solutions. 
ATLAS uses the R\RAMESH (iMacNeice et al.ll200Qh hierarchi- 
cal block-structured adaptive mesh refinement for high eff'ective 
resolution in areas of physical interest. The solenoidal constraint 
(V • B = 0) is preserved using a staggered mesh algorithm based 
on the lBalsara & Spiced (il999) fiel d transport method. ATLAS 
uses a Piecewise Parabolic scheme dColella & Woodw^l 19841) 
to reconstruct the values at cell inter faces and the Roe-Balsara 
app roximate MHD Riemann solver (lRodll98i1: iRoe & Balsaral 
[1996) to explicitly compute the cell interfaces fluxes. The multi- 
dimensional correction used to compute the transverse fluxes is 
the Corner Transport Upwind scheme of lColellal (Il99d) as mod- 



Domain 

Refinement 
Max Resolution 
N Jet deprojected velocity 
S Jet deprojected velocity 
Jet radius 
Jet density 
Ambient density 
Jet Temperature 
Ambient Temp 
Angle between jets 
Angle between jets 



-840AU<3;,z<840AU 

0<x<1680AU 

5 levels 

2.625 AU 

Vj = 300 km s-i 

Vj = 100 km s-i 

rj = lO^^cm 

rij = 500 cm-3 

ria = 5000 cm-3 

Tj = lO'^K 

Ta = lOOTT 

(Case I & III) 

10 (Case II) 



Table 3. Initial conditions for simulation 



ified by lSaltznianl (Il994l) . In our own case, the simulations were 
carried out using ATLAS on a 64 node cluster of the so-called 
"Beowulf" type. 

3.4. Initial and boundary conditions 

Consistently with the location of the object in Taurus, we as- 
sumed a distance of 140 parsecs to the jets. We used values 
derived from the observations for the density, temperature and 
velocities, as described in Table 2. The ambient medium is mod- 
elled with a uniform density (^« = 5 x 10^ cm"^) and tempera- 
ture (10^^) and the jets are modelled with density rij = O.lria, 
tempe rature lO'^K a nd velocities of 100 and 300 km s"^ respec- 
tively (iLiseau et al. 20 05). Regarding t he flow velocity, we take 
the radial velocities of iHartiganet al.l (l2QQQb) 200 km s"^ and 
60 kr n s~^ and deproiect them from the plane of the sky assuming 
(^with lHartigan et alT(l2000b l)) an inclination angle of 45 degrees. 
In order to reproduce the knotty structure of the flows we also 
use a sinusoidally varying injection velocity (Raga et al. 1990), 
assuming an amplitude of ±30% in the velocity with a period of 
8 years for each jet. We also stagger the launching of the two 
jets, arguing that the faster northern jet is launched 150 years af- 
ter the slower southern jet. The calculations of the launch times 
are based on estimating the age of the jets from their lengths 
and current velocities, assuming the velocities have not changed 
over time. The velocity profile across the jet section is a posi- 
tive cosine - with its maximum at Vjet and its minimum (zero) at 
Yjet. This profile is chosen based on the high angular re solution 
observations of the DG Tau jet by lBacciotti et al.l (12000 ). where 
the highest velocities appear to be found in the axial region of 
the jet. 

Typical postshock c ooling lengths estimat ed using the plane 
parallel shock models of JHartigan et al.l (Il987l) are of the order of 
30 AU in the colliding outer wings of the binary jet bow shocks 
and in the internal shocks so will be resolved in ~12 points. 

The initial jets' diameters are approximately 13 points each 
and the jet cross-sectional areas are resolved in approximately 
150 square cells for each jet. 

The boundary conditions are inflow for r < r^, reflecting for 
r > rj along the inner x-boundary (x=0) and outflow along all 
other boundaries. 



4, Results 

We present here results from fully 3D HD and MHD simulations 
of the jets. 
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4.1 . Case I: Hydrodynamic case 

By using a simple hydrodynamical (HD) model, we were able to 
reproduce the twisted morphology of the LI 551 IRS 5 outflow, 
that in this case appears to arise from the interaction of the two 
jets. In the observations of L1551 IRS 5, the two jets emerge 
from a binary source separated by 45 AU. The observed angular 
separation is 20 degrees which corresponds to a deprojected an- 
gle of approximately 10 degrees. Initially, for the hydrodynamic 
case we assume the jets are parallel. (In the next section (14.21) we 
remove this simplifying assumption.) We model the jet propaga- 
tion and show that the interaction afl'ects the propagation. 



300 
250 
200 
150 
100 

50 


50 100 150 200 250 300 



i'^^^I^^^HHH^I^^^HI^B^I^^^^^H 



-0.050.38 0,80 1.23 1.66 2,08 2.51 



Fig. 3. Midplane slice colour density map of 3-D binary jet sim- 
ulation at ^ = 190 years. The scale of the grid is 1680AU by 
1680AU. 




1.769 2.006 2.245 2.480 2.717 2.954 5.191 




500 



1000 



1500 2000 



921 2.151 2.582 2.612 2.842 5.075 




500 1000 1500 2000 



1.291 1.624 1.957 2.290 2.625 2.956 5.289 




500 1000 1500 2000 



0.815 1.241 1.668 2.096 2.524 2.951 5.579 



The density slice from the 3D HD simulation (see Figure [3]) 
clearly shows that the slow, secondary jet appears to bend close 
to the inlet. Such a bend or kink is also observed in the slow, 
southern jet of L1551 IRS 5 about 4'' (56 AU at a distance 
of 140 pc from the source (lltoh et al.ll2000l) ). In the simulation 
shown in Figure [3] a kink in the slow southern jet is also visible 
- both in the density midplane cut shown in Figure [3] and in the 
derived emission maps shown in [71 This is caused by the bow 
shock of the fast jet interacting with the beam of the slow jet. 
The northern jet has a Mach number ~ 3 times higher than the 
slow southern jet and simply pushes it out of the way. This re- 
produces the observed kink at 4'' - without the need for magnetic 
fields. There is no noticeable reaction by the fast jet - possibly 
indicating that the estimated velocity is too high. 

4.2. Case II: Magneto-hydrodynamic case 

In the second case, we consider the inclusion of the magnetic 
field and run a full magneto-hydrodynamic (MHD) simulation, 
trying to reproduce the jet behaviour closer to the source. The 
magnetic field provides a mechanism for the southern jet to 
change its direction - hence we allow the jets to initially diverge. 



Fig. 4. Time evolution of log jet total density in the MHD bi- 
nary jets at 1=42.5,85,127.5,170 years. The images shown are 
2D midplane cuts through the 3D grid. 



The projected angle betw een the two jets is about 20 degrees 
(tFridlund & Liseaul Il99 8) we use a deprojected angle of 10°. 
To redirect the two jets would require either a density contrast 
e.g. th e wall of a conven iently shaped cavity or a magnetic hoop 
stress te allv et al.ir2003h . 

We us e the magnetic field configuration based on the obser- 
vations of lScarroti(ll988h . that suggest the presence of a toroidal 
magnetic field around the system. The toroidal field may be pro- 
duced by the circumbinary disk twisting the froz en-in local mag- 
netic field lines (see, e.g. IMatt & Pudritzl 120041 and the refer- 
ences therein.) 

We assume that the ambient medium is permeated by a 
toroidal field with a maximum value of B= 10 yuG. The magnetic 
vector potential A, where B = Vx A, had the analytical form A^^ 
cos 271 X -\- cos 271 y. This is an approximation to the real field 
which would be helical according to the disk-twisting theory. 
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.769 2.006 2.243 2.480 2.717 2.954 3.191 




200 400 600 800 



.691 1.921 2.151 2.382 2.612 2.842 3.073 




400 600 



800 



1.701 1.960 2.219 2.478 2.737 2.996 3.255 




400 600 



800 



0.993 1.313 1.634 1.954 2.274 2.595 2.915 



Fig. 5. Section of Figured] showing velocity vectors. 



lower ram pressure and thus is easier to redirect. In the simu- 
lation over the time evolution we see the southern jet change 
its direction to move parallel to the axis of the toroidal field. 
Hence, it appears that the hoop stress from the jet toroidal field 
slowly changes the direction of the jets. Thus, magnetised binary 
jets will tend to collimate and refocus along the direction of the 
fastest or strongest jet. 

4.3. Case III: Orbiting binary jet 
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Fig. 6. Surface plot of density of 3-D binary jet simulation at 
t = 5, 10, 40, 60, 90 years. The orbital period of the source is 260 
years. The jet pair has performed 1/3 of a rotation. 



Figure |4] shows a series of density slices from the 3D MHD 
simulation. The binary outflow is modelled in full 3D with 
toroidal magnetic field. Compared to the previous HD simula- 
tion, we have set the angle between the two jets to be close to 
10°. In the hydrodynamic case there is no mechanism to to force 
the jets to become parallel. In the MHD case the Lorentz force 
can redirect both of the jets. As a consequence, we would expect 
to see the two jets propagating into two diff'erent directions on 
large scales. On the contrary, the smaller jet is refocused along 
the faster jet. The secondary jet (i.e., the south low velocity one), 
bends towards the z-axis, impelled by Lorentz force. This can be 
seen very clearly in FigureHl panel 1 for instance, where the slow 
jet is already moving towards the centre before the approach of 
the fast jet. Eventually, the jet material from the slower jet gets 
merged into the faster one. 

The observations of LI 551 show both jets diverging from 
the source until about 4'' from the source when they change di- 
rection to pursue a roughly parallel course. The change in di- 
rection is most pronounced for the southern jet which has the 



We now explore the eff'ects of orbital motion of the jet 
sources on the survival of the binary jets. iMasciadri & Ragal 
(2002) have modelled jets from orbiting sources However they 
modelled just a single jet orbiting a central body and compared 
its evolution to an analytical prediction. Figure [S] shows a model 
of a binary jet from two sources orbiting a common centre of 
mass. The jet crosses the grid in 90 years. In this case, for sim- 
plicity we assume that the jets are launched at the same time and 
their velocities are 300 and 100 km s"^ 

We find that, by calculating the dynamical age of the jets in 
LI 551 IRS 5, the orbital period of the binary -260 years is too 
long to have much eff'ect on the fast northern jet (crossing time 
-90 years). 

4.4. Emission maps 

In order to properly compare the results of our simulations with 
the existing observations we have computed synthetic emis- 
sion maps in the light of the forbidden lines of neutral oxygen 
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([OI]/16300). In order to predict the line emission produced by 
the two jets in our model we use the density, temperature and 
the fraction of ionised hydrogen as computed by the numerical 
simulation, to find the emissivity in each point of the jet gas. The 
emissivity e is defined as the energy emitted at the wavelength A 
of the transition by a unit volume of gas per unit time. Following 
iBacciottil (l2QQ2h , the emissivity can be expressed as: 



— '^-^dKDfe) 



(7) 



where A^ is the transition probability, hh is the hydrogen nu- 
clear number density, Xe is the ionised fraction of hydrogen, ^ 
is the fraction of the element Z ionised at level i, ^ is the abun- 
dance of element Z with respect to hydrogen, and ^^ is the 
fractional population of the upper level of the transition. We as- 
sume s olar abundances for the elern ents in the fluid. For O, fol- 
lowing |Bacciottr^Eisi5^1 (|1999|) . we can determine the frac- 
tional ionisation as a function of Xg assuming the relationship 



O^ 



(Co + So)Xe 



O^ (ao - So)xe + S'^ 



(8) 



holds, where Co is the collisional ionization rate, ao the direct 
plus dielectronic recombination rates. So and S'^ are the direct 
and inverse charge ionisation exchange rates respectively. The 
fractional population of the upper level, ^^^ , is basically a func- 
tion of the free electron density and the electron temperature, and 
is determined by a 5 -level linear system of equations describing 
the statistical equilibrium of the level populations subject to col- 
lisional excitation and de-excitation, and to spontaneous radia- 
tive dec ay. A more cornplete d escription of the met hod can be 
found in 'Bacci otti et al.l(ll995l) : lBacciotti & Eislo"feil(fl999.) and 
iBacciot ti (2002). 

The emissivity values calculated in the 3D space from our 
simulation are then integrated along lines of sight, assuming an 
inclination angle of 45 degrees. The result is then convolved with 
the beam of the instrument, which corresponds to an angular res- 
olution of 0.'' 1 in the HST case. 

The synthetic emission map produced in this way is shown in 
Figure[7l It is clear that the interaction has a strong observational 
signature. The second jet is virtually obliterated but remains vis- 
ible at the point where the fast jet's bow shock impinges on its 
beam. The emission map is quite difl'erent to what is seen in the 
density slice. The knotty structure of the southern jet and the 
asymmetric structure of the northern jet are only apparent in the 
emission map. These may be compared with the observations in 
Figure [H Even though the image has been convolved to account 
for the beam resolution a peak in emission is visible in the south- 
ern jet where the bow shock of the fast jet is striking the beam. 
In the observations a similar peak in the emission is visible near 
the source of the southern jet (Knot Z in Figure [T]). 



5. Discussion 

The LI 551 IRS 5 jet pair can be seen as a pathological object 
unlike any other observed at such close range. In fact many ef- 
fects, like binary source orbiting, peculiar environment, density 
parameters, magnetic fields play a role in sculpting this most 
unusual object. To attempt an explanation of the observed phe- 
nomena that can take into account all the above mentioned ele- 
ments, a full 3D numerical treatment is necessary. In this paper 
we have presented 3D HD and MHD simulations of interact- 
ing binary jets that may help explain the L 1551-IRS5 case. We 
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Fig. 7. Synthetic [OI]/16300 emission map derived from the 
physical parameters of our simulation in the hydrodynamic case 
as depicted in Figure[3l to be compared with corresponding HST 
observations of the L1551 IRS 5 jet (HH154) (See text) The im- 
age shows the northern and southern jets at a time t=190 years. 
Knot Z indicates the point where the bow shock of the northern 
jet collides with the beam of the southern jet. The units are erg 
cm"^ s"^ arcsec"^. 



have used the code ATLAS including radiative cooling and non- 
equilibrium ionisation treatment. The main results can be sum- 
marised as follows: 

- We can reproduce the kink and bending of the secondary jet 
and the knotty structure of the two components as observed 
forL1551IRS5. 

- Over a long time-scale, we see precession induced by the 
orbital motion of the source. On the short lifetime of the jets 
from L1551, this eff'ect is negligible. 

- If the jets are not strictly parallel, as in most observed cases, 
we show that the magnetic field can help the collimation and 
refocusing of both the jets. In fact the toroidal field aff'ects 
the motion of the southern jet; and this magnetically driven 
change in direction together with the interaction of the bow 
shock of the fast jet with the beam of the slow jet contribute 
to the distinctive morphology. 

- We have produced emission maps in the [01] lines which 
can be compared with observations. We show that the kink 
structure of interacting jets is still apparent in the synthetic 
observations. 

Some aspects however remain to be clarified. For example, 
one of them is the nature of the the source of the X-ray emission. 
Bally et al. (2003) find a source of X-rays in L1551 IRS 5 which 
they attribute to either fast shocks or reflected x-rays scattered 
out through the outflow cavity. The fast shocks may be caused by 
jet collimation which could be magnetic in nature. In our model 
there is a peak in emission (Knot Z in Figure [B which appears 
to come from the colliding winds. The Bally model suggests that 
the X-rays come from a moving source at the base of the jets. If 
the density ratio is riambient I ^ jet =10 then the shock speed will 
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be too slow to achieve X-ray emission temperatures. If on the 
other hand the density ratio is lower i.e. 0.1 the jet is moving 
into a less dense region swept out by the slower jet the shock 
velocities will be higher. To get up to the observed shock velocity 
the density contrast would need to be much lower than observed 
e.g. 0.01 would give a velocity in the range required. We also 
could assume that it is caused by magnetic reconnection in the 
interval between the two jets where the field is compressed by 
the pair of bow shocks. 

One other interesting aspect to be discussed is the implica- 
tions for jet launching of the binarity of the source. The binary 
provides a vastly diff'erent natal environment for a protostar and 
its attendant jets and outflows than that envisaged in a single star 
formation. We identify three main difl'erences for jet launching 
for binaries, orbital motion, circumbinary material and interac- 
tion with a jet. In this work we showed that the impact upon the 
jet beam of the bow shock of its neighbour is large enough to 
be observable as a knot of emission moving along along the jet. 
The orbital motion is of too long a period to afl'ect the jet prop- 
agation. The role of the circumbinary material may be to twist 
up the magnetic field into the toroidal shape which appears in 
the polarized light. In conclusion, our simulations have demon- 
strated to be able to reproduce the observed signatures in this 
source and to give powerful constraints on the physical eff'ects at 
play. A larger number of observed cases has however to be mod- 
eled with our code in order to give a complete view of binary jet 
production. 
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